Non-Hermitian anisotropic XY model with intrinsic rotation-time reversal symmetry 
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We systematically study the non-Hermitian version of the one-dimensional anisotropic XY model, 
which in its original form, is a unique exactly solvable quantum spin model for understanding 
the quantum phase transition. The distinguishing features of this model are that it has full real 
spectrum if all the eigenvectors are intrinsic rotation-time reversal (7?.T) symmetric rather than 
parity-time reversal ('PT)symmetric, and that its Hermitian counterpart is shown approximately to 
be an experimentally accessible system, an isotropic XY spin chain with nearest neighbor coupling. 
Based on the exact solution, exceptional points which separated the unbroken and broken symmetry 
regions are obtained and lie on a hyperbola in the thermodynamic limit. It provides a nice paradigm 
to elucidate the complex quantum mechanics theory for a quantum spin system. 
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I. INTRODUCTION 



In recent years, much effort has been devoted to es- 
tabhsh a parity-time (VT) symmetric quantum theory 
as a complex extension of the conventional quantum 
mechanicsi^ since the seminal discovery by Bender—. 
A cornerstone of the theory is the fact that the non- 
Hermitian Hamiltonian with VT symmetry can have an 
entirely real energy spectrum. It is also motivated by the 
interest in complex potentials in both theoretical and ex- 
perimental aspects, since the imaginary potential could 
be realized by complex index in optics^"^^. A natural 
question to ask is whether the non-Hermitian quantum 
theory prefers to the VT symmetry. While there is as 
yet no answer to this question, we can gain some insight 
regarding the pseudo-Hermiticity of a non-Hermitian 
model without the VT symmetry. On the other hand, 
the non-Hermitian quantum model in discrete system is 
a nice testing ground for the study of the non-Hermitian 
quantum mechanics because of its analytical and numer- 
ical tractability. As simplified models, tight-binding and 
quantum spin models can capture the essential features of 
many discrete systems. In general, the non-Hermiticity 
arises from VT symmetric on-site imaginary potentials 
in tight-binding models^ '^^i^^" — , and imaginary magnetic 
fields in quantum spin model o^^'^^ . In addition, the com- 
plex coupling constant between spins can also introduce 
the non-Hermiticity, which relates to the spin, the in- 
trinsic degree of freedom. Then the spin rotation may 
take the role of the parity operation for constructing a 
pseudo-Hermitian spin system. 

In this paper, we propose a pseudo-Hermitian model 
without VT symmetry explicitly, but with intrinsic 
rotation-time (TZT) reversal symmetry. It is a non- 
Hermitian version of the one-dimensional anisotropic XY 
model. The original Hermitian XY model was initially 
solved and has became the paradigm of the quantum 
spin system possessing a second-order quantum phase 
transition^*. We will show that this model has full real 
spectrum if all the eigenvectors have TZT rather than 
parity-time reversal (VT) symmetry, and its Hermitian 



counterpart is an isotropic XY spin chain. The result for 
such a concrete example may have profound theoretical 
and methodological implications. 

This paper is organized as follows. In Section |TT1 we 
present the model Hamiltonian. Based on the solutions 
we investigate the phase diagram and analyze the sym- 
metry of the ground state. In Section Hill we construct 
the Hermitian counterpart of the model and its approx- 
imate reduced form. Finally, we give a summary and 
discussion in Section HVl 



II. 



HAMILTONIAN AND INTRINSIC TZT 
SYMMETRY 



The model Hamiltonian of the non-Hermitian 
anisotropic one-dimensional spin-i XY model in a trans- 
verse magnetic field A for N particles is given by 



N 



H 



y y , \ z 



(1) 



y, z) are the Pauli operators on site j, and 



satisfy the periodic boundary condition ct" = cr"+Af- -'^'^'^ 
the sake of simplicity, we only concern the case of even N , 
the conclusion is available in the case of odd TV. In com- 
parison with the model proposed by Giorgi^^, the present 
model is not VT symmetric. The non-Hermiticity arises 
from the imaginary anisotropic parameters ±17. The 
Hermitian version of the XY model is completely solved 
by applying the Jordan- Wigner— , Fourier and Bogoli- 
ubov transformations2M2. The Jordan- Wigner transfor- 
mation maps the Pauli operators into canonical fermions, 
while the Fourier transformation essentially decomposes 
the Hamiltonian into invariant subspaces due to the 
translational symmetry of the system. We will see 
that these transformations are applicable in solving the 
Hamiltonian ([T]) if we extend it to its complex versions. 
Before solving the Hamiltonian, it is profitable to inves- 
tigate the symmetry of the system and its breaking in 



2 



the eigenstates. By direct derivation, we have [TZ, H] 
and [T,H] 7^ 0, but 



to replace the Pauli operators by the fermionic operators 
Cj . We note that the parity of the number of fermions 



[TZT, H]^0. 



(2) 



The Hamiltonian is rotation-time (Ti-T) reversal invari- 
ant, where the linear rotation operator TZ has the function 
of rotating each spin by 7r/2 about the z-axis 



TZ — cxp 



(3) 



and the antilinear time reversal operator T has the func- 
tion TiT = —i. Before we consider the general non- 
Hermitian XY model, we highlight key ideas on two lim- 
iting cases. When 7 = 0, the Hamiltonian ((TJ reduces to 
the ordinary XY model with external field in z direction 



N 



(4) 



which has full real spectrum and its eigenstates can al- 
ways be written as the eigenstates of operator TZT- On 
the other hand, when 7 ^ A and 1, the Hamiltonian ([1]) 
reduces to 



N 



(5) 



which has full imaginary spectrum. Then for an eigen- 
state \ipn) of Hoc with nonzero eigenvalue En, i.e., 

Hao \i>n) = En l-^n) , we have 



HooTZT\i'n 



^EnTZTl^n 



(6) 



due to the facts of Eq. © and = The eigen- 

state \tpn) obviously breaks the TZT symmetry. These re- 
sults strongly imply that the TZT symmetry in the present 
model plays the same role as T'T symmetry in the T'T 
pseudo-Hermitian system. It motivates further study of 
such a model systematically. 



N 



(8) 



is a conservative quantity, i.e., [i?, H] — 0, where Np 



^f^ic'jCj. Then the Hamiltonian ([T} can be rewritten 



as 



where 



P, = -(l-|-r;H) 



(9) 



(10) 



is the projector on the subspaces with even (77 = +) and 
odd {n — —) Np. The Hamiltonian in each invariant 
subspaces has the form 



-'7 (cjyCi -I- cjcjv -I- i7cjyci -I- i7CiC7v) (11) 

N 

-2JX^c]cj+NJ\. 
Taking the Fourier transformation 



-y 



(12) 



for the Hamiltonians we have 



A. Solutions 

Now we consider the solution of the non-Hermitian XY 
Hamiltonian of Eq. ([1]). We note that the general solu- 
tions of the Hermitian anisotropic XY Hamiltonian are 
not restricted to the real anisotropic parameter. We will 
show that in the case of imaginary parameter, eigenstates 
and energies are still accessible. As the same procedures 
performed in solving the Hermitian Hamiltonian, we take 
the Jordan- Wigner transformation^ 



a| = l-2c]c„ 



(7) 



-7sinfc„ (c_fc,Cfc^ + cLfc^c[,J - A] 



(13) 



where k+ = 2tt {m + 1/2) /N , fc_ = 2TTm/N, m = 
0,1,2,...,7V-1. 

So far the procedures are the same as those for solving 
the Hermitian version o{ H. To diagonalize the non- 
Hermitian Hamiltonian, we introduce the Bogoliubov 
transformation in complex version: 

= cos (?) Ck„ - i sin (^j cLfc,^ (14) 



Ak^ = cos ( ^ 1 



2, 

2 ) ""k^ 



2 sni I - I c-k„ 
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where 



tan (0) = 



Z7 sin 
(A — cos krj) 



(15) 



We would like to point out that this is the crucial step to 
solve the non-Hermitian Hamiltonian, which essentially 
establish the biorthogonal bases. Obviously, complex Bo- 
goliubov modes (Afe,, , A^^ ) satisfy the canonical commu- 
tation relations 

{A,,,A,,} = 4,„fc;,, (16) 

which result in the diagonal form of the Hamiltonian 

fi„=^e(A,fc„7)(^A^A^-i^. (17) 

Here the single-particle spectrum in each subspace is 

e (A, kjj, 7) = — 2jy^ (A — cosfc^)^ — 7^ sin^ k^j. (18) 

Note that the Hamiltonian is still non-Hermitian due 
to the fact that Ak^ ^ A^ . Accordingly, the eigenstates 



of H,^ can be written as the form 



IG 



')/ J 



(19) 



which constructs the biorthogonal set associated with the 
eigenstates 



of the Hamiltonian ^ where 



n 

fc>0 



(20) 



|Vac) (21) 



is the ground state of iJ^, and |Vac) is the vacuum state 
of the fermion Cj. In the following, we will investigate the 
phase diagram based on the properties of the solutions. 



B. Phase diagram 

It is clear that when any one of the momentum A;,, 
satisfies 



|A — cosfc^l < |7sinfc,, I 



(22) 



the imaginary energy level appears in single-particle spec- 
trum, which leads to the occurrence of complex energy 
level for the Hamiltonian ([1]), and the TZT symmetry is 
broken in the corresponding eigenstates. This can be 
seen from the properties of the single-particle spectrum 
and the ground states |G,,). 

Firstly, we focus on the boundary between the broken 
and unbroken symmetry regions. There are totally 2N 
equations in the form of |A — cosfc,, | = |7sinA;,, | for all 
the possible value of fc,,, only one of which determines 
the line along the boundary of the diagram within certain 
region in the parameters A and 7 plane. Then the phase 
boundary is dog-leg path for finite N, and becomes a 
smooth loop for the infinite N. In the thermodynamic 
limit, the momentum fc,, becomes continuous, then the 
phase boundary as a curve can be given by the following 
set of parametric equations^^ 



0, 



de (Ac,fc^,7c) 
dkjf 

e (Ac, fc,,,7c) = 0. 



(23) 
(24) 



Straightforward algebra gives the analytical boundary 
curve as 



Ic 



1, 



(25) 



which is a hyperbola. In addition, the broken region 
does not include the line 7 = 0. In Fig. [1] we plot 
the phase diagrams for the systems with = 4, 8, and 
30, respectively. It is shown that as N increases, the 
boundary approaches to the hyperbola of Eq. (1^51) . 

Secondly, according to the non-Hermitian quantum 
theory, the occurence of the exceptional point always ac- 
complishes the TZT symmetry breaking of an eigenstate. 
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FIG. 2. (Color online) Plots of the coupling constant (a) k(0), (b) k(1), (c) k (2) and (d) k (3) for the system of N = 10. The 
blue solid line (empty triangle), red solid line (empty circle) and black solid line (empty square) indicates the approximately 
analytical p9|) (numerical ((35])) result for 7 = 0.2, 7 = 0.5 and 7 = 2, respectively. It is shown that the approximately analytical 
results are in agreement with the numerical results as A increases. Note that the scales of the subfigures for various coupling 
constants k are different. 



For the concerned model, the symmetry of the ground- 
state |G^) can be indicator of the phase transition due 
to the fact that the groundstate energy becomes complex 
once the system is in broken region. In the following, we 
focus on the discussion about the symmetry of |G^) in 
the different regions. 

Taking the combination of the Jordan- Wigner and 
Fourier transformations on the rotational operator in Eq. 
([3]), we have 



(26) 



r;=±,fc^ 



where = cj, Ck^ is the particle number in fc^ space. 
Applying the TZT operator on the fermion operators and 
its vacuum state |Vac), we have 



and 



TZTcl (7^r)-^ = ^ct_,^, 



Ten Vac) = |Vac) , 



(27) 



which are available in the both regions. However, the 
coefficients cos(6'/2) and sin (0/2) experience a transi- 
tion as following when the corresponding single-particle 
level changes from real to imaginary: We have [cos (6/2)]* 
= cos(6'/2) and [sin(6'/2)]* = -sin(6'/2) for real levels 
and [cos (0/2)]* = sin (0/2) for the imaginary levels, re- 
spectively. This leads to the conclusion that the ground- 
state is not TZT symmetric in the broken region, i.e., 



nT\Gr,) = \Gr,) ; Unbroken 
7^r|G,,) 7^ |G,,) , Broken 



(29) 



(28) 



It shows that the TZT symmetry in the present model 
plays the same role as VT symmetry in the VT pseudo- 
Hermitian system. 

As a comparison, it is noted that the phase boundary 
of the Hermitian anisotropic XY model is an ellipse. It 
is worth pointing out that the phase transitions in the 
Hermitian and the non-Hermitian models are the dif- 
ferent types of quantum phase transition. The typical 
quantum phase transition^i describes an abrupt change 
in the ground state of a many-body system. For the Her- 
mitian anisotropic XY model, the transition occurs when 
the order parameter, the energy gap between the ground 
and first excited states, goes to zero. For the present 
non-Hermitian model, zero gap also leads to the phase 
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boundary of the ground state, which tends to — 7^ = 1 
in thermodynamic hmit. It accords with the boundary of 
the unbroken TZT symmetric region. It reveals two differ- 
ences between the phase transitions in the non-Hermitian 
and the Hermitian XY modeh (i) The energy gap can 
vanish for the former with the finite N , while for the 
Hermitian XY model zero-gap is never achieved unless 
the thermodynamic limit is reached, (ii) There is no TZT 
symmetry breaking in the quantum phase transition of 
the Hermitian XY model. 



III. HERMITIAN COUNTERPART 

In the complex quantum theory, as an important 
perspective, the physical meaning of a non-Hermitian 
Hamiltonian has received a lot of attentions. When 
speaking of the physical significance of a non-Hermitian 
Hamiltonian, one of the ways is to seek its Hermitian 
counterparts^, which possess the identical real spec- 
trum. According to the complex quantum mechanics, 
a non-Hcrmitian Hamiltonian can be transformed into 
a Hermitian Hamiltonian by introducing a metric, a 
bounded positive-definite Hermitian operator—. However, 
the obtained equivalent Hermitian Hamiltonian is usu- 
ally quite complicated^ii^. On the one hand, it is tough 
to provide an explicit mapping of a pseudo-Hermitian 
Hamiltonian to its equivalent Hermitian Hamiltonian 
and on the other hand, the obtained Hermitian coun- 
terpart is unphysical, which usually contains the nonlo- 
cal interactions^ii^. In this section, we will provide a 
physical counterpart of the present model and analyze 
its behavior when approaches to the exceptional point. 

A Hermitian Hamiltonian becomes a Hermitian coun- 
terpart of a non-Hermitian system if and only if they 
share the identical real spectrum. In this sense, the coun- 
terpart is not unique and seems to be easy constructed. 
However, it is expected that the obtained counterpart is 
physically relevant, or possessing the local interaction, 
which makes things difficult. To this end, our strategy is 
to consider the operator as the form 



We employ the Fourier transformation 



(30) 



where dk^ is Fermion operator satisfying 

{d,^,d,,} = {4^,4,}=0. 



(31) 



Here /i+ and /i_ represent the Hamiltonians in the in- 
variant subspaces with even number of '^1+^'=+ ^^'^ 
odd number of ^ J, 4 , respectively. Obviously, /i^ 
is Hermitian when e (A, fc,p 7) is real and has the identical 
spectrum with that of i/,,. 



dk. 



and Jordan- Wigner transformation^^ 

Kj 



(32) 



(33) 



to rewrite the Hamiltonian by the Pauli spin operators 
t" {a = X, y, z). Here {t"} and {erf} are different sets 
of spin operators, which have no direct relation, or more 
precisely, there is no explicit transformation to connect 
them. The essence is that we replace the pair of conjugate 

operators ^^fc^ , A^'^ ^ by (^d^. , dk'^ ^ , which ensures the 
identical spectrum of two Hamiltonians. However, they 
are distinct due to the fact that {Ak^ ^ Ak^ ■ 

Now we transform the Hamiltonians hn into the fol- 
lowing spin model 

h„ = -^Y.Ki^ n ^/^m^r+H.c. + ^KnTf, (34) 



l<j<m 



where 



J 



(I — m) 



(35) 



is T^-dependent coupling constant. Here /i+ and rep- 
resent the Hamiltonians in the invariant subspaces with 
even and odd number of magnons. It looks compli- 
cated due to the long-range coupling and the extra phase 
ni<j<m ^'^^^ show that it can be reduced approx- 

imately to a simple model in the non-trivial parameter 
region. 

First of all, we consider the strong filed limit case, A ^ 
y^l + 7^. Taking the Taylor expansion, we get 



k^(|Z-to|) « — ^[(A-cosfc^) 

(A + cos kn) 7^ sin^ ^ ^,(\i-m\)k„ 
2A2 

Using the relation 
we have 

3 

Kjj {\l - m|) = [k (n) (5|/_„i|,„ + O (i/^)] 

ra=0 



(36) 



(37) 



(38) 
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>.=2, 7=0.1 




FIG. 3. (Color online) Plots of the energy spectrum for the 
systems of = 4, 6, 8 with A = 2, 7 = 0.1. The black and 
azure solid lines denotes the energy spectrum of Hamiltonian 
H and T-i, respectively. It can be observed that the energy 
spectra obtained by two Hamiltonians are in well agreement 
with each other. 



feature of k,, {I — m) as (A, 7) turns to the critical point 
('^olc) except the case of Ac = 1 (7c = 0), which is not 
the exceptional point. To this end, we take the gradient 
of K,j {I — m) in the A — 7 plane, i.e.. 



dX 



ei 



97 



■62 



(41) 



^ Y: cos [k, [l - m)] " '''' ^-^''-^ ^"^^ 



\J {\ — cos fc)))^ — 7^^ sin^ kr, 



where ei and 62 denote the unit vectors, • Cj = 5ij. 
When (A, 7) turns to the critical point (Ac, 7c), the 
dominant contribution to the summation is the term of 
kri — kc, where 



I. 1 

COSKc = — . 

Ac 



(42) 



Then we have 



where f = 7^/A^ and 



k{0) 



1 



AJ, 



2 V 8 
«(2) = ^,«(3) 



(39) 



It is shown that {I — m) is N and 7y-independent, and 
decays as |? — to| increases. Secondly, numerical simu- 
lation indicates that even the field A is not so strong, 
k{1 — m) is still negligible for |Z — m| > 1. Fig. [5] is the 
plot of K (Z - m) in Eq. ((35)) . 

Then the Hamiltonian can be reduced approximately 

as 



lim Vk 

(A,7)->(A,,7,) 



(Ac — cos k„) ei — 7c sin fc„e2 
X ' ' - = 00. 

y (Ac - cos fc^)^ - 72 sin^ fc,, 

It accords with the fact that the derivative of the energy 
diverges when the system tends to the exceptional point. 
Then for a Hermitian matrix, the non-analytic behavior 
is not from the Jordan block but from the divergence of 
derivatives of matrix elements. The similar situation also 
occurs in another modeP^, which may imply a general 
conclusion. 



IV. SUMMARY 



n^-^Yl i^^^r+i + H.c.) + (0) ^ rf , (40) 
I I 

which is the ordinary XY model. To demonstrate this 
equivalent, we make a comparison of H and the original 
Hamiltonian H in Eq. (jlj. The energy levels are com- 
puted for small size systems. In Fig. [3] we see that they 
are in agreement well in the physical parameter region. 

Secondly, it is worthy to point out that H does not con- 
tain the information of the exceptional point if we take 
the approximate form of k (0) and k (1) in Eq. As 
well known, distinguished from phase transition in a Her- 
mitian system, the non-analytic properties for the phase 
transition of a non-Hermitian are caused by the Hamilto- 
nian becoming a Jordan block operator at the exceptional 
point. Nevertheless, the Jordan block cannot appear in 
a Hermitian Hamiltonian. Therefore it is interesting to 
see what happens to the equivalent Hamiltonian at the 
exceptional point. To investigate the critical behaviors 
of the equivalent Hamiltonian, the exact coupling con- 
stants {I — m) should be considered. We focus on the 



In summary, we have proposed a non-Hermitian model 
without VT symmetry explicitly, but with intrinsic TZT 
symmetry, which is a non-Hermitian version of the 
anisotropic XY model. Based on the exact solution, we 
have found the phase diagram for the finite N system and 
the corresponding symmetry breaking in ground state. 
An analysis of the symmetry in the ground state shows 
that the TZT symmetry is broken when the groundstate 
energy becomes complex, which is similar to that for a 
VT system. It indicates that the VT symmetric model 
is not the unique candidate of pseudo-Hermitian system. 
We also constructed a Hermitian counterpart which ex- 
hibiting the identical real spectrum to that of the orig- 
inal one within the unbroken region. It is an isotropic 
XY spin chain with long-range coupling, which has been 
shown to have the following properties: The Hermitian 
counterpart can be reduced to a simple XY model with 
nearest neighbor coupling when the system is not in the 
vicinity of the exceptional point. The derivatives of all 
the coupling strengths with respect to the system pa- 
rameters diverge at the critical points. This provides an 
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interpretation for the critical behavior of level repulsion. 2012CB921900. 
The result for such a concrete example may have pro- 
found theoretical and methodological implications. 
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